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Summary 

This paper investigates the application of the linear-quadratic- Gaussian 
(LQG)/loop transfer recovery (LTR) method to the problem of synthesizing an 
attitude control system for a large flexible space antenna. The study is based 
on a finite-element model of the 122-m hoop/column antenna, which consists 
of three rigid-body rotational modes and the first 10 elastic modes. A robust 
compensator design for achieving the required performance bandwidth in the 
presence of modeling uncertainties is obtained by using the LQG/LTR method 
for loop shaping in the frequency domain. Based on the limitations of achieved 
performance with the LQG/LTR method, different sensor actuator locations 
are analyzed in terms of the pole zero locations of the multivariable systems and 
possible best locations are indicated. The computations are performed by using 
the LQG design package ORACLS augmented with frequency domain singular 
value analysis software. 


1. Introduction 

The robust control synthesis methodology which 
uses frequency domain matrix norm bounds (i.e., sin- 
gular values) has received considerable attention in 
the recent literature (refs. 1 through 5). 1 The ba- 
sic framework for frequency domain synthesis using 
the linear-quadratic-Gaussian (LQG) /loop-transfer- 
recovery (LTR) methodology was developed in refer- 
ences 1 and 2. It has been applied to diverse systems 
such as power systems (ref. 3) and aircraft engine 
control (ref. 4). The LQG/LTR design philosophy 
uses a low frequency “design model” of the plant 
and a high frequency characterization of the mod- 
eling errors. This method, which uses unstructured 
uncertainty and singular value bounds, appears to 
be particularly well suited for the control of large 
flexible spacecraft because of the considerable un- 
certainty that inherently exists in the mathematical 
models. 

Large space structures (LSS) are typically highly 
flexible and can be represented by partial differen- 
tial equations or by very large systems of ordinary 
differential equations. They have many resonant fre- 
quencies, some of which may be very low and possi- 
bly closely spaced. The natural damping is usually 
very small. For these reasons, control of LSS is a 
challenging task (ref. 5). 

Since the system is inherently of high order, a 
practical controller has to be based on a reduced- 
order design model. Furthermore, the parameters 
(i.e., frequencies, mode shapes, and damping ratios) 

1 This methodology was also presented in a lecture entitled 
“Multivariable Control System Design Using the LQG/LTR 
Methodology” given by Professor Michael Athans of the 
Massachusetts Institute of Technology at the Langley 
Research Center in 1984. 


of the system are known imprecisely; this introduces 
additional modeling errors. Typically, the modeling 
errors for finite-element models increase substantially 
with increasing modal frequency. Controller design 
for large space structures using frequency domain 
analysis was investigated in reference 6, wherein sin- 
gular value analysis was used to check the robustness 
with respect to modeling errors. 


One of the planned activities of the NASA Space 
Transportation System is the placement in Earth or- 
bit of a variety of large space antennas. Potential 
large space missions will require antennas and struc- 
tures ranging from 30 m to 20 km in size. Ap- 
plications include communications (mobile, trunk- 
ing, etc.), remote sensing (soil moisture, salinity, 
etc.), deep space network (orbital relays), astronomy 
(X-ray, observatory, optical array, radio telescope, 
very long baseline interferometry, etc.), energy, and 
space platforms. Specific missions have been pin- 
pointed and future requirements identified for large 
space antennas for communications, Earth sensing, 
and radio astronomy (ref. 7). Particular emphasis 
is placed on mesh-deployable antennas in the 50- to 
120-m-diameter category. One such antenna is the 
maypole (hoop/column) antenna, shown schemati- 
cally in figure 1. The hoop/column antenna basi- 
cally consists of a central mast attached to a de- 
ployable hoop by cables in tension (refs. 7 and 8). 
The deployable mast consists of a number of tele- 
scoping sections, and the hoop consists of 48 rigid 
segments. The reflective mesh, which is made of 
knit gold-plated molybdenum wire, is attached to the 
hoop by graphite fibers. The mesh is shaped by us- 
ing a network of stringers and ties to form the radio- 
frequency (RF) reflective surface. In order to achieve 



required RF performance, the antenna must be con- 
trolled to specified precision in attitude and shape. 
For example, for missions such as land mobile satel- 
lite system (LMSS) for providing mobile telephone 
service to users in the continental United States, it is 
necessary to achieve a pointing accuracy of 0.03° rms 
(root mean square) and a surface accuracy of 6 mm 
rms. It is also necessary to have stringent control 
(usually a fraction of a degree) on the motion of the 
feed (located near one end of the mast) relative to 
the mesh. Because of its large size and relatively 
light weight, the antenna is highly flexible, having a 
large number of significant elastic modes. It is there- 
fore necessary to use reduced-order controllers. 

Reduced-order control synthesis for this system 
using the standard LQG theory was investigated in 
reference 9. The standard LQG procedure yielded 
satisfactory control, that is, rigid-body bandwidth of 
up to 0.25 rad/sec, satisfactory time constants for the 
elastic modes, and acceptable rms pointing errors in 
the presence of sensor noise. It should be noted that 
the LQG approach in reference 10 used a large num- 
ber of actuators and sensors (four three-axis torque 
actuators and four three-axis attitude and rate sen- 
sors). It was found that the first three flexible modes 
had to be included in the design model (in addition 
to the three rigid modes) to obtain satisfactory per- 
formance. The main problem with the LQG method 
was that a large number of weighting parameters had 
to be simultaneously adjusted to obtain a good de- 
sign. In addition, the stability robustness proper- 
ties with respect to inaccuracies in the modal pa- 
rameters could not be properly assessed because it 
was difficult to effectively characterize the bounds 
on modeling errors in a time-domain setting. The 
LQG/LTR method offers an approach for overcom- 
ing these problems in the frequency domain. 

The purpose of this paper is to investigate the 
use of LQG/LTR multivariable frequency domain 
methodology in the design of an attitude control sys- 
tem for the hoop/column antenna. The design ap- 
proach considered herein is determiiiistic, and the 
rms pointing errors in the presence of noise are not 
analyzed. A sequence of low order compensators is 
obtained by considering design models which include 
increasingly higher frequency modes. The final de- 
sign model chosen is the first one in this sequence 
which allows the performance/ robustness objectives 
to be met. In this sequence of design models, the 
first one consists of the rigid-body modes only. Sub- 
sequent design models are obtained by the successive 
addition of flexible modes. The designs use three- 
axis torque actuators, collocated attitude and rate 
sensors, and attitude and rate feedback. The organi- 
zation of this paper is as follows. The mathematical 


model of the system is first described in section 3. 
Section 4 describes the control objectives, followed 
by a brief description of the LQG/LTR technique 
(refs. 1 and 2) in section 5. The reduced-order (low 
frequency) design model and the high frequency mod- 
eling uncertainty barrier are also discussed in that 
section. The numerical results are given in section 6. 
Section 6.1 presents the results of synthesizing the 
controller based on the LQG/LTR procedure with 
only attitude feedback being used. Section 6.2 de- 
scribes the results obtained when both attitude and 
rate feedback are used for controller design. Some of 
the problems and limitations observed are also high- 
lighted. Specifically the locations of transmission ze- 
ros and their impact on the performance for two dif- 
ferent sensor and actuator locations are discussed. 
The results are summarized in section 7. 

2. Symbols 

A system matrix 

A p system matrix for full-order plant 

B input matrix 

B p input matrix for full-order plant 

C output matrix 

C p output matrix for full-order plant 

D damping matrix 

G regulator gain matrix 

G c compensator transfer function 

G c modified compensator function 

G f complete plant transfer function 

G^ Kalman filter transfer function for 
attitude feedback 

Grf Kalman filter transfer function for 
attitude and rate feedback 

G^ modified loop gain 

G p plant transfer function 

G p modified plant transfer function 

Gi rigid-body transfer function 

G 2 elastic-mode transfer function 

H Kalman gain 

I identity matrix 

I k k X k identity matrix 

l s inertia matrix 

L weighting matrix used in Kalman filter 

design 
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L p multiplicative uncertainty 

l m upper bound on multiplicative 

uncertainty 

m order of input vector 

n order of system 

n q number of elastic modes 

rtx number of torque actuators 

P regulator Riccati matrix 

q modal amplitude vector 

q regulator design parameter 

r reference input 

s Laplace variable 

T j torque vector for j th actuator 

u input vector 

x state vector 

x state estimate vector 

y a sensed attitude 

y r sensed rate 

a s rigid-body attitude vector 

AG additive uncertainty 

& damping ratio for ith elastic mode 

9 S rigid-body rotation about Y-axis 

A diagonal matrix of squares of natural 

frequencies 

/i Kalman filter design parameter 

£ Kalman filter Riccati matrix 

a singular value 

a largest singular value 

a smallest singular value 

mode slope matrix at actuator 
locations 

4> s rigid-body rotation about X-axis 

mode slope matrix at sensor locations 
i>s rigid-body rotation about Z- axis 

(jj frequency 

natural frequency of ith elastic mode 

Dots over symbol indicate derivative with respect 
to time; superscript T indicates a transposed matrix. 


3. Mathematical Model 

As a consequence of its large size and light weight, 
the antenna is a highly flexible system having a large 
number of significant structural modes. A finite- 
element model of the antenna (ref. 8) is used in this 
paper. The mathematical model considered consists 
of rotational rigid-body dynamics (about the three 
axes) and the elastic motion. We assume that the 
control will be accomplished by using nx three-axis 
torque actuators. The linearized equations of motion 
are 


nj> 

Is«s = T i 
3 = 1 

(1) 

q + Dq + Aq = 

(2) 


where I s is the 3x3 inertia matrix, T j is the three- 
axis torque applied by the j th actuator, 
ot s — {(frsiOsi'tps) denotes the rigid-body attitude, 
q is the n q x 1 modal amplitude vector (for n q struc- 
tural modes), D = 2 diag (ftwi, £2^2, • . , $nq<*>nq) 
is the inherent damping matrix (where & is the damp- 
ing ratio for the ith mode), <& is the 
m x n q “mode-slope” matrix (where m = 3 nj>), 
u = (T|\ . . . , T^ t ) t is the m x 1 vector of 

actuator torques, and A — diag (wj, . . . , oj^ q ) 
where is the frequency of the ith elastic mode. 
The rigid-body parameters and the first 10 elastic 
frequencies are given in table I. Each value of Q is 
assumed to be 0.01 for i = 1, 2, . . . , n q . 

Normally, the sensors used include attitude and 
rate sensors. The three-axis attitude y a at a sensor 
(e.g., a star tracker) location is given by 

ya = 0L S + (3) 

where 9 is the 3 X n q mode-slope matrix at the sen- 
sor location. If an attitude rate sensor (e.g,, a rate 
gyro) is used, the sensor output y r is given by an 
equation similar to equation (3), except that a s and 
q are replaced by 6t s and q, respectively. The torque 
actuators and the attitude and rate sensors are as- 
sumed to be located near the top and bottom of the 
mast as shown in figure 1. The elastic deformation 
caused by some of the flexible modes is shown in fig- 
ure 2. 

Define x = q^\q T ) ; then the state 

space model can be written in the following form: 

x = A^x + B^u (4) 

y a = C F x (5) 

The sensor and actuator noise are not considered 
in the design process in this paper; however, it 
will have to be included when computing the rms 
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pointing errors. Ignoring the noise, the transfer 
matrix between the input (three-axis torque) and the 
output (three-axis attitude) is given by 

Gf( s ) = Gi(s) + G 2 (s) (6) 


assumes that the loop is broken at the plant output. 
Adjust the weighting matrices in the KBF design un- 
til its frequency response meets the robustness spec- 
ifications at high frequencies and bandwidth specifi- 
cation at low frequencies. 


where 


I -1 

Gi(s) = 

s 


n q 

° 2 ^ § S 2 + 2 + u. 




(7) 

(8) 


where and represent the mode-slope matrices 
at the sensor and actuator locations corresponding 
to the zth mode. Numerical data used to define 
equations (4) through (8) and a discussion of the 
model characteristics are given in the appendix. 


4. Design Objectives 

The basic design objectives for the control sys- 
tems are (1) to obtain sufficiently high bandwidth 
(i.e., closed-loop frequencies corresponding to the 
rigid-body modes) and satisfactory closed-loop damp- 
ing ratios for rigid-body and structural modes; and 
(2) to obtain satisfactory rms pointing errors, feed 
motion errors, and surface errors. The first design 
objective arises from the need to obtain sufficiently 
fast error decay when a step disturbance (such as sud- 
den thermal distortion caused by entering or leaving 
the Earth’s shadow) occurs. The second design ob- 
jective arises from the RF performance requirements. 
These two objectives may not necessarily be com- 
patible. For example, the use of increased feedback 
gains for obtaining higher bandwidth and damping 
ratios will in general, result in higher rms errors (be- 
cause of the amplified effect of sensor noise) beyond 
a certain point. Therefore, it is necessary to carefully 
consider the trade-offs between the speed of response 
and lower rms error. In this study, the main control 
system specification is to ensure a minimum band- 
width of 0.1 rad/sec for the closed-loop system. The 
upper limit on the low frequency gain is not specified 
but it is desired that it should be as high as possible. 


5. Design Procedure 

The LQG/LTR method (refs. 1 through 3) basi- 
cally consists of the following steps: 

1. Define a “design” model of the nominal plant 
which is an acceptable low frequency representation. 
Define the high frequency uncertainty (robustness) 
barrier and the low frequency performance barrier. 

2. Design a full state feedback compensator based 
on the steady-state Kalman-Bucy filter (KBF). This 


3. Design an LQ regulator to asymptotically “re- 
cover” the frequency response obtained in step 2. 

4. Verify stability, robustness, and performance 
for the entire closed-loop system. 

The first step, which consists of the definition of 
the plant and the uncertainty (robustness) barrier, 
is often the most important one. The basic prob- 
lem in controlling a flexible structure is the pres- 
ence of a large number of lightly damped struc- 
tural modes. Practical limitations necessitate the 
use of reduced-order controllers. Therefore, the un- 
controlled modes and the error in the knowledge of 
the controlled modes together represent uncertainty. 
Since the number of structural modes is usually large 
and finite-element modeling accuracy typically de- 
creases with increasing model frequency, the design 
model should consist of rigid-body plus the first few 
elastic modes. In order to obtain an acceptable low 
frequency representation, the design model must in- 
clude at least the three rigid-body modes. The uncer- 
tainty barrier is a measure of the plant uncertainty 
at high frequencies. The plant uncertainty can be 
represented as either multiplicative or additive un- 
certainty (fig. 3). The LQG/LTR approach requires 
the characterization of the uncertainty in terms of 
a frequency-dependent upper bound. Frequency do- 
main sufficient conditions are used to test the ro- 
bustness in the presence of uncertainties within that 
bound. 

For the case of multiplicative uncertainty, L p (joj) 
of figure 3(a), with the loop broken at the plant 
output, the closed-loop system is stable if 

a [Lp(jw) - I] < o jl + [Gp(j'w) G c (jw)] _1 | (9) 

where G p (joj) and G c (jcu) are the design model 
(plant) and compensator transfer matrices, and a 
and a denote the largest and the smallest sin- 
gular values of the argument matrix, respectively. 
At high frequencies, assuming ||Lp(ju/)|| 1 and 

l|Gp(ju>) G c (jc^)|| 1 (ref. 8) (approximately) 

yields 

<t(G p G c ) < l/cr(Lp) (10) 

•For notational convenience, the argument u joj n has 
been dropped. 

The uncertainty (or robustness) barrier is an up- 
per bound lm(w) on a(L p ). The system is stable 
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in the presence of such unstructured uncertainties if 
<7[G P G C ] < l m 1 ( oj) at high frequencies. 

When the additive uncertainty formulation 
(fig. 3(b)) is used, a sufficient condition for stabil- 
ity robustness is given by (ref. 10) 


cr(AG) < 


a(I + G P G C ) 

^(Gc) 


(U) 


At high frequencies, assuming ||G P G C || <C 1 (ref. 10) 
(approximately) yields 


^(Gc) < 


1 

cr(AG) 


( 12 ) 


That is, the compensator must roll off sufficiently 
rapidly at high frequencies. The main feature of 
the LQG/LTR approach is to first design a full state 
compensator (based on KBF) which has the behavior 
of the desired loop transfer matrix (i.e., the loop gain 
GpG c ). Therefore (from step 2), any loop shaping 
should involve the product G p G e rather than G c 
alone as in equations (11) and (12). Assuming that 
G p is a square matrix, 


G c — Gp 1 (G P G C ) (13) 

^(Gc) < <r(Gp x ) a(G p G c ) 


The advantage of an LQG-based full state design is 
that it has excellent classical properties, and its fre- 
quency response can be shaped in the desired manner 
by varying the weighting matrices (ref. 2). As dis- 
cussed in reference 2, this design can be accomplished 
by using the linear-quadratic regulator (LQR) Ric- 
cati equation if the loop is broken at the plant in- 
put, or the KBF Riccati equation if it is broken at 
the point where the residual signal enters the KBF. 
Herein we select the latter because the objective is 
to control the attitude output. This selection is also 
consistent with references 3 and 4. The KBF equa- 
tions are 


AS + SA t + LL t - SCrCS = 0 


(16) 




(17) 


where L and // are the design parameters, L being 
an n X m matrix and fj, a scalar, H is the KBF gain 
matrix, and S is the corresponding Riccati matrix. 
The KBF loop transfer matrix is given by 


G/ff(s) = C(sl - A) -1 H (18) 


or 

<7(G C ) <£ -1 (G p )cf(GpG c ) (14) 

Using equations (14) and (11), the following sufficient 
condition for stability robustness is obtained: 


<t(I + GpG c ) a(Gp) 
cr(GpG c ) 


> <r(AG) 


(15) 


Generally, the frequency response <t[GxfU^)] would 
shift higher as fx decreases, and the crossover fre- 
quency can be adjusted by changing L. 2 

Having obtained satisfactory singular value be- 
havior of KBF, the third step is to design an LQR to 
“recover” the desired frequency response. This is ac- 
complished by solving the algebraic Riccati equation 


The second step in the design procedure is to de- 
sign a full state feedback compensator having desir- 
able singular value properties. The performance of 
the closed-loop system depends on the low frequency 
gain and the crossover frequency of the loop transfer 
matrix GpG c , that is, on the behavior of a[GpG c ]. 
Larger low frequency gain and crossover frequency in- 
dicates better tracking performance. Thus, a[GpG c ] 
should lie above the performance specification as 
shown in figure 4(a). The other requirement is the 
stability robustness in the presence of model uncer- 
tainties. If the multiplicative uncertainty formulation 
is used, according to equation (10), the a[GpG c ] plot 
should pass under the robustness barrier a~ l (L p ) at 
high frequencies (fig. 4(a)). On the other hand, if 
the additive formulation is used, the robustness con- 
dition, equation (15), should be satisfied (fig. 4(b)). 


A t P + PA - PBB r P + qC T C = 0 (19) 

where P is the Riccati matrix and q is a positive 
scalar. It has been proven in references 1 and 2 
that the loop transfer matrix G P G C for the overall 
system (consisting of the plant, the KBF, and the 
LQR) tends to G^ j p(s) as q — > oo, provided that 
the open-loop plant has no transmission zeros in the 
right half-plane. After recovery, G c (s) is given by 

G c (s) = G(sl - A + BG + HC) _1 H 


2 Ibid. 
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A block diagram of the compensator is shown in 
figure 5. 

Since the compensation obtained has no guaran- 
teed robustness properties, the last step will consist 
of testing the eigenvalues of the entire closed-loop 
system to ensure the stability and robustness. The 
overall closed-loop system is 


x 

x 


Ap 

-b f g 

x 

nc F 

A-BG-HC 

X 


where the subscript F is used to denote the full-order 
nominal plant, and x denotes the state estimate for 
the design model. If instability is discovered, it will 
be necessary to return to step 2 and redesign the KBF 
for lower bandwidth and the LQR for robustness 
recovery. If this does not produce satisfactory results, 
it would then be necessary to return to step 1 and 
include more elastic modes in the design model. 
Application of the foregoing LQG/LTR procedure 
for the hoop/column antenna is described in the 
following section. 

6. Numerical Results 

Before applying the LQG/LTR method for com- 
pensator design, the controllability, observability, 
and invariant zero properties of the plant (with re- 
spect to different sensor/actuator locations) were 
studied. Based on these studies (presented in de- 
tail in the appendix), location 1 shown in figure 1 
was selected for the sensors and the actuators. The 
LQG/LTR procedure was subsequently applied, first 
by using only attitude sensors and then by using at- 
titude and attitude-rate sensors. The computations 
were performed by using the LQG design package 
ORACLS (ref. 11) augmented with frequency domain 
singular value analysis software. 

6.1. Controller Design Using Attitude 
Feedback 

In this section, the results of the application of 
the LQG/LTR procedure are presented when only 
attitude feedback is employed. The studies are 
performed by using the design model consisting of 
(1) only rigid-body model (n = 6, n q = 0), (2) rigid- 
body and the first flexible mode (n = 8, n q — 1), 
and (3) rigid-body and the first three flexible modes 
(n = 12, n q =3). The measurements available are 
the three attitude angles at location 1. One three- 
axis torque actuator is used at the same location. 
The compensator is designed based on these sensors 
and actuators. The results of the designs are pre- 
sented in figures 6 through 23. 

The largest and the smallest singular values of 
the rigid-body transfer matrix (n = 6) are plotted 


in figure 6. The corresponding additive uncertainty 
AG, which consists of the (20th order) flexible dy- 
namics, is plotted in figure 7. Figure 6 clearly shows 
the roll-off behavior of the plant, which is of the 
order of 1/s 2 . Figure 7 indicates the presence of 
poles near the undamped flexible mode frequencies 
of 0.75 rad/sec, 1.35 rad/sec, etc. Also, the pole near 
the first mode frequency near 0.75 rad/sec produces 
the highest peak. (The importance of this fact is seen 
later when the stability condition is violated at this 
point.) It can also be seen that a zero is present near 
5 rad /sec. 

Based on this rigid design model, the following 
choice of L and // were made after a number of trials: 


L = 


O3 

io- 2 i 3 


fl = 1 


The objective is to obtain a bandwidth of at least 
0.1 rad/sec. The resulting values of the “desired” 
loop gain matrix G are shown in figure 8. 

The standard LQG/LTR procedure requires the 
definition of the desired loop transfer characteristics 
(see step 2 in section 5). That is, a( G rp) must 
satisfy the low frequency performance specifications, 
and <f( Gkf) must satisfy the high frequency robust- 
ness specifications. Thus, in the presence of additive 
uncertainty AG, the procedure states that the ro- 
bustness condition 

a{I + Gxp)a(G p ) . 

— " ff(AG) 

should be satisfied. However, in the computations 
performed, it was found that this condition makes 
the desired design (G kjp) extremely conservative in 
the sense that higher loop gains can be used without 
causing instability. Therefore, recovering this con- 
servative loop gain yields a compensator with poor 
performance. This fact led to a modification of the 
LQG/LTR procedure. In particular, this robustness 
test on Gkf is omitted in the modified procedure. 
Instead, the recovery is carried out first, and then the 
(less conservative) stability test (eq. (11)) is applied 
directly for the compensator G c . 

With the use of the recovery procedure, the com- 
pensator is obtained for this case with q — 10 4 . The 
resulting stability test (eq. (11)) is shown in figure 9. 
It is seen that the stability margin is lowest at the 
first mode frequency (0.75 rad/sec). Any increase in 
the gain (obtained by q > 10 4 ) resulted in the vio- 
lation of the stability condition at that point. The 
resulting singular values of the compensator G c are 
given in figure 10, which shows that the compen- 
sator is of the “lead-lag” type. The overall loop 
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bandwidth is obtained from the singular values of the 
loop transfer function G p G c shown in figure 11. It is 
seen that the bandwidth (i.e., the frequency at which 
o;(GpG c ) — 1) is far short of the required 0.1 rad/sec. 
The plots roll off as 1/s 2 (40 dB/decade) and indi- 
cate the recovery process. In order to increase the 
bandwidth, the regulator gain has to be increased by 
increasing q . However, this results in the violation 
of the stability condition (eq. (11)). Thus it became 
evident that, with the rigid-body design model, it is 
not possible to meet the performance specifications. 

To overcome these problems, the inclusion of the 
first flexible mode (0.75 rad/sec) in the design model 
was considered. The inclusion of the first flexible 
mode, which is predominantly a torsion mode, results 
in a design model of order 8. The singular value 
plots for G p and AG are given in figures 12 and 
13. Figure 12 indicates the presence of a pole of G p 
near 0.75 rad/sec and a zero near 0.082 rad/sec. The 
singular values of AG in figure 13 are an order of 
magnitude lower than those in figure 7 (wherein AG 
consisted of all the flexible modes). The first peak of 
a(AG) occurs at 1.35 rad/sec, which is the frequency 
of the second mode. This is the critical point in the 
stability test (eq. (11)). After a number of trials, the 
following choice of L and fi was made: 

O3 

Kr 2 i 3 

io~ 2 i 2 j 0 

The singular values of Gpp are plotted in fig- 
ure 14. The recovery is obtained for q = IQ 5 and the 
stability test is shown in figure 15. 

Figure 15 indicates the critical point to be about 
0.28 rad/sec. There is a good margin at the peaks of 
AG due to upward sloping of the upper curve. The 
resulting compensator is shown in figure 16. The 
compensator has a pole at about 0.25 rad/sec and 
a zero at 0.8 rad/sec. The resulting loop transfer 
function (G p G c ) plots are shown in figure 17. The 
plots indicate that the required 0.1 rad/sec band- 
width is not obtained (although it is much higher 
than the rigid-model case). Any increase in the gain 
(for q > 10 5 ) was found to result in the violation 
of the stability condition (eq. (11)). Figure 17 in- 
dicates the presence of the open-loop invariant zero 
near 0.082 rad/sec, which was also confirmed by in- 
dependent computations. This zero is almost on the 
imaginary axis (i.e., the transfer matrix is close to 
being nonminimum phase). Therefore (as would be 




expected), the recovery procedure is not very effec- 
tive for making GpG c approximate Gj^p. 

In order to improve the performance further, the 
next step was to include the first three flexible modes 
in the design model. It is logical to do this because 
they represent the first modes about each axis, that 
is, the first torsion mode and the first bending modes 
in the X-Z and Y-Z planes. Thus the order of 
the design model was 12. The singular value plots 
for Gp and AG are shown in figures 18 and 19, 
respectively. Figure 18 indicates that Gp has zeros 
near 0.082 and 0.22 rad/sec and poles near 0.75, 1.35, 
and 1.7 rad/sec. It is seen from the AG plot (fig. 19) 
that <x is considerably lower than that in figures 7 and 
13. After numerous trials, the following choice of the 
L matrix and the scalar fi were made: 

O3 

10-U3 


10" 4 I 2 I 0 
10- 4 I 2 I 0 
10- 4 I 2 I 0 

The resulting G^ with the desired characteris- 
tics is shown in figure 20. The recovery was accom- 
plished with q = 10 10 . The stability test is shown in 
figure 21. It can be seen that condition (11) is satis- 
fied with a wide margin. Also, at the peaks for AG 
(at 8 rad/sec), the upper curve slopes upward, indi- 
cating good tolerance of high frequency uncertainty. 
The limit for increasing the gain (indicated by the 
lowest point in the upper curve in figure 21) occurs 
at about 0.3 rad/sec. The resulting compensator G c 
is shown in figure 22. The gain of G c is much higher 
than that obtained in the previous cases. Generally, 
the LQG/LTR technique attempts to choose G c in 
such a way that Gp is approximately inverted and the 
product GpG c is replaced by Gxp. The three- mode 
design plant shown in figure 18 has elastic mode 
eigenvalues at -0.0075 ± y0.75, -0.0135 ± j 1.35, 
and —0.017 ± jl.7. Figure 22 shows that G c has ze- 
ros with frequencies near these locations. The design 
model also has transmission zeros at 

(-0.9 x 1(T 4 ) ± yo.082 
(-0.37 x 10 -3 ) ± j0.22 
(-0.29 x 10 -3 ) ± j0.22 
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The zeros of the design model are too close to 
the juj axis and tend to numerically behave as 
nonminimum phase. Some alleviation is obtained 
by the compensator pole near 0.4 rad/sec. The 
plots for the loop transfer matrix GpG c are given 
in figure 23. It is seen that a bandwidth of 
0.1 rad/sec is obtained except for the presence 
of the invariant zero near 0.082 rad/sec which 
causes some deterioration of performance. This 
zero is invariant under feedback and depends on the 
sensor/actuator locations. At higher frequencies, 
G P G C rolls off at 60 dB/decade. Also, a and q_ are 
closely spaced; this indicates good system behavior. 
The compensator gains (i.e., the regulator and filter 
gain matrices) for this case are shown in table II, and 
the complete closed-loop eigenvalues are shown in 
table III. 

Thus it is seen that the inclusion of the first three 
modes in the design model yields a robust compen- 
sator which also meets the bandwidth specifications. 


For the design with attitude and rate feedback, 
the transfer matrix for the design model is 


G p (s) 


’ G S s ) ■ 

.sG(s). 


6x3 


where G(s) is the transfer function matrix with 
attitude as the output. Thus G p (s) is always of 
rank 3, and v(G p G c ) = 0. Therefore the loop gain 
can never satisfy the bandwidth specifications. To 
overcome this difficulty, we consider only the attitude 
measurements for performance evaluation. In other 
words, the compensator G c uses both the attitude 
and the rate as its inputs (and has the dimension 
3x6), but we consider only the 3x1 attitude in the 
performance evaluation. Figure 24(b) shows the loop 
with a slight modification to account for this. The 
loop gain for the system is 

G l (s) = [I 3 0 3 ] [G P G C ] [l 3 sh} T 


6.2. Controller Design Using Attitude and 

Rate Feedback 

In this section, the numerical results for the de- 
sign using the feedback of both attitude and rate are 
discussed. Similar to the attitude feedback case, one 
three-axis torque actuator is used at location 1. One 
three-axis attitude sensor and one three-axis attitude 
rate sensor are assumed to be located at the same lo- 
cation (i.e., location 1). Thus the plant G p (s) is now 
6x3, with the three attitude angles and the three 
attitude rates being the outputs and three torques 
being the inputs. The compensator G c (s) is then 
3x6. The resulting loop is shown in figure 24(a). 
The realization of G c consists of a KBF using the 
six measurements, followed by the 3 x n LQ regu- 
lator gain matrix G, where n is the order of the de- 
sign model. As in the attitude feedback design, three 
cases are considered, with the design model consist- 
ing of 

Case a: Only rigid-body modes (n = 6) 

Case b: Rigid-body and the first flexible mode 
(n = 8) 

Case c: Rigid-body and the first three flexible 
modes (n = 12) 


and one needs to check ^(G^) as the measure of 
performance (i.e., bandwidth). Similarly, in the KBF 
design for the desired loop gain, one needs to check 
a of 

Gkf = [I 3 O 3 ] G kf [!3Sl3] T 
As in the previous case, the robustness test (eq. (11)) 
is carried out by using (the 6x6) GpG c (and not 
G kf) although a(G^) is used as the measure of 
performance. 

The design procedure involves first solving the 
KBF equations with suitable choices of L and fi 
to obtain the desired performance characteristics 
^(Gkf)* The next step is to recover this desir- 
able loop gain by increasing q the LQ regulator de- 
sign. However, the loop transfer recovery can be ob- 
tained only if the number of inputs is greater than 
or equal to the number of outputs, which is not the 
case here. Based on the experience with attitude 
feedback designs, q cannot be usually increased suf- 
ficiently to obtain recovery in the desired frequency 
range (because the stability condition is violated). 
Therefore, although asymptotic recovery is not possi- 
ble, a decent bandwidth should still be obtainable by 
increasing q . 

Figure 25 shows the frequency response of G p for 
case a, i.e., rigid-body design model. The a - plots of 
AG are shown in figure 26. After a number of trials, 
the following L and were used: 


As in the design with attitude feedback, the second 
and third modes (being bending modes about X- and 
F-axes) were not considered separately, and hence a 
three-mode model was considered after the one-mode 
model. 


03 

03 

1 

L °3 ! 

1 

CO 

1— 1 

1 

O 

tH 


/i = 1 


J 6x6 
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The resulting cr-plots of Gkf are shown in figure 27. 
After the KBF design, the recovery procedure is 
carried out by increasing q . The largest q without 
causing instability (as indicated by the stability test 
(eq. (11)), which is shown in fig. 28) was q = 10 4 . 
Figures 29 and 30, respectively, show the cr-plots of 

Gl and G c = G c [I35I3] . It is seen from figure (29) 
that the bandwidth requirement of 0.1 rad/sec is 
practically impossible to meet using a rigid model. 

For case b, the design model consisted of the rigid- 
body and the. first flexible modes. Figures 31 and 32 
show the a-plots for G p and AG, and figure 33 shows 
the plots for Gkf which was obtained with 

O 3 j O 3 

o 3 | 10-U3 


10 X l2 | 02x4 




in figure 40. The a-plots of Gi and G c are shown 
in figures 41 and 42. As shown in figure 41, the 
bandwidth requirement is not satisfactorily met. In 
addition, the presence of a transmission zero near 
0.082 rad/sec introduces a notch in the cr(Gjr,) plots 
similar to the attitude feedback case. The com- 
pensator gains and the closed-loop eigenvalues of 
the complete system are given in tables IV and V, 
respectively. 

Based on these results, it appears that designing 
with attitude feedback only allows for higher perfor- 
mance bandwidth while maintaining stability robust- 
ness. The reason can be seen from a comparison of 
the stability robustness condition (eq. (11)) in both 
the attitude and attitude and rate cases. In the atti- 
tude case, AG 3 X 3 rolls off as 1/s 2 . For attitude and 
rate, it rolls off as 1/s. Since, in both situations, G c 
behaves as 1/s and a(I+G p G c ) — ► I for large values 
of s, condition (11) is ultimately satisfied. However, 
comparing figures 21 and 40 shows that the attitude 
and rate case experiences greater difficulty in the fre- 
quency region containing the peak of <7 (AG). It can 
be shown that 


The largest q in the LQ regulator design (which 
does not cause instability) was q = 10 5 , and the 
corresponding stability test (eq. (11)) is shown in 
figure 34. The a-plots of G/, and G c are shown 
in figures 35 and 36, respectively. It is clear that 
the bandwidth specification cannot be met with this 
model. 

In order to improve the performance, the first 
three flexible modes were next included in the model 
(in addition to the rigid modes). Figures 37 and 38, 
respectively, show the cr-plots of G p and AG. The 

Gkf plot shown in figure 39 was obtained for the 
following choice of L and (x\ 


O 3 j O 3 

o 3 | 10 - : % 

lO-Ua | 0 2 x4 ** = 1 

10 U 2 j 02x4 

, I 

10 l l 2 i 0ox4 

z I j 12x6 

The largest permissible q (without causing instabil- 
ity) was 10 8 , and the resulting stability test is shown 



^(AG6x3(j^)) — (1 + a; 2 ) 1 / 2 <7(AG3 X 30 W )) 

The (1 + u; 2 ) 1 / 2 magnifies a(G(j(jj)) and causes a 
lower value of q in the recovery process and a more 
conservative performance bandwidth. This difficulty 
could possibly be overcome by a different choice of L 
and /i. 

7. Concluding Remarks 

The LQG/LTR (linear-quadratic-Gaussian/loop 
transfer recovery) multivariable frequency domain 
technique was employed in the design of an attitude 
control system for large flexible space antenna. The 
LQG/LTR method was noted to be especially attrac- 
tive for overcoming spillover effects common to large 
space structures control problem modeled by using 
finite-element data. The design objective of avoiding 
excitation of higher order modes while satisfying per- 
formance criteria was met by including these modes 
in the robustness uncertainty barrier. 

Design was based on a reduced-order model cho- 
sen as the rigid-body dynamics plus the fewest num- 
ber of low frequency vibrational modes necessary to 
meet a desired closed-loop bandwidth. Inclusion of 
the first three vibrational modes (corresponding to 
the three axes) was found to be necessary to meet 
a bandwidth of 0.1 rad/sec. For higher bandwidths, 
design models with greater than three modes would 
be needed. A satisfactory control design was ob- 
tained by using only a single collocated pair of three- 
axis attitude sensor and torque actuator for this 
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problem. An investigation employing attitude-rate 
sensors in addition to the attitude sensors showed no 
improvement over this arrangement. 

Performance degradation was observed due to 
the presence of invariant zeros within the design 
bandwidth. These zeros were unavoidable given the 
prescribed sensor/actuator locations and emphasized 
the fact that consideration should be given to control 
aspects right from the early design phases of large 
space structures. 


A modification of the standard LQG/LTR pro- 
cedure was introduced in which the robustness test 
was performed with the the full LQG compensator in- 
stead of the intermediate Kalman filter design. This 
approach was found to produce higher gain compen- 
sators and helped overcome the basic conservative- 
ness shortcoming of the LQG/LTR approach. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
April 15, 1986 
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Appendix 

System Characteristics of Hoop/Column 
Antenna Mathematical Models 

In this appendix, the state-variable realizations 
employed in the LQG/LTR design studies are ana- 
lyzed to determine their innate controllability, ob- 
servability, and invariant zero properties. Compu- 
tational results are presented for a single three-axis 
torque actuator and a collocated attitude sensor. A 
combined attitude and rate sensor case yielded the 
same results. 

The hoop/column antenna model given by equa- 
tions (4) and (5) are 

x = A^x + B^u (Al) 

y a = C F x (A2) 

In equation (Al), with n — 6 + 2 n q) 


Equation (A4) describes the plant rigid-body 
motion, and equation (A5) describes the elastic 
motion contribution when n q vibrational modes 
are included in the design model. Frequencies 
Ui(i = 1, 2, ..., 10) are given in table I for a 
122-m-diameter hoop/column structure. The damp- 
ing ratios £ are taken as & = 0.01(z = 1, 2, . . . , 10). 
Eigenvalues of A z E (i == 1, 2, . . . , 10) are as given in 
table AL Employing a single three-axis torque actu- 



where, in units of (ft-Ib-sec 2 ) 1 , 

I^ 1 = diag (2.37 x 1CT 7 , 2.37 x 1CT 7 , 3.09 x 10 -7 ) 

(A8) 

The matrix V g depends on actuator location. 
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For actuator location 1 (fig. 1): 


Vb = 


For actuator location 2 (fig. 1): 




0 

6.3 x 10“ 5 
0 
0 
0 

3.71 x 1(T 3 
0 

-7.49 x 10“ 7 
0 

-6.21 x 10“ 4 
0 

5.55 x 10“ 6 
0 
0 
0 

1.1 x 10“ 3 
0 
0 
0 

2.43 x 10“ 3 


0 

-1.81 x 10 -5 
0 
0 
0 

-3.76 x 1(T 4 
0 

1.70 x 10- 7 
0 

-1.08 x 10" 2 
0 

6.69 x 1(P 6 
0 
0 
0 

4.41 x 1(P 3 
0 
0 
0 

_ 4.90 x 1(P 3 


0 

0 

0 

-2.91 x 1(P 3 
0 
0 
0 
0 
0 
0 
0 
0 
0 

-5.00 x 1CP 4 
0 
0 
0 

-3.04 x HP 3 
0 
0 


0 

0 

0 

2.52 x 1(P 4 
0 
0 
0 
0 
0 
0 
0 
0 
0 

-6.47 x 10" 3 
0 
0 
0 

-7.58 x 10“ 3 
0 
0 


0 

-5.04 x 10“ 3 
0 
0 
0 

5.45 x 1(P 4 
0 

-8.45 x HP 6 
0 

-5.6 x 1(P 5 
0 

-3.19 x 1(P 4 
0 
0 
0 

-2.23 x 10~ 4 
0 
0 
0 

-1.05 x 10“ 4 


0 

-1.62 x 1(P 5 
0 
0 
0 

1.22 x 10“ 5 
0 

-1.97 x 1(P 4 
0 

-4.09 x 10~ 6 
0 

-1.26 x 10“ 2 
0 
0 
0 

-1.23 x 10~ 5 
0 
0 
0 

-2.16 x 10“ 6 


}mode 1 
}mode 2 
}mode 3 


}mode 10 


}mode 1 
}mode 2 
}mode 3 


}mode 10 


(A9) 


(A10) 


In both equations (A9) and (A 10), mode slope mag- 
nitudes less than 1CP 9 have been set to zero. 

For collocated attitude sensors and torque 
actuators, 

Cf = [(Is, 0), Vdsxn ( An ) 

with 

Vc = [ UVslLlh, (A12) 


U< = 


0 

1 



(A14) 


The matrix in equation (A12) is chosen from 
equation (A9) or (A10), depending on actuator loca- 
tion. For collocated rate sensors, 

Cf-[(0,I»),VS] Jxn (A15) 


U 2nqx2nq — diag (Uj, U2 } * * • » U n ^) (A13) 


A measure of the relative controllability of the 
rigid-body and elastic modes can be obtained through 
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examination of the reciprocal condition numbers of 
the polynomial matrix (ref. 12) 

C(s) = [(sI n -A F ),B F } (A16) 


over the set of eigenvalues of A^. Reciprocal con- 
dition numbers are defined as the ratio of smallest 
to largest singular values for C(s). Comparison of 
these numbers evaluated at s = 0.0 for the rigid- 
body mode and one of each of the eigenvalues of 
A^, {i — 1 , 2 ,..., n q ) gives an indication of the rel- 
ative ease or difficulty in affecting a particular mode 
through the control input. A similar approach ap- 
plied to the polynomial matrix 


0(s ) 



(A17) 


can be used to measure relative observability prop- 
erties. Table All presents controllability and observ- 
ability reciprocal condition numbers for collocated 
sensor and actuator locations 1 and 2 in figure 1. 

Table All indicates that of all the modes con- 
sidered, the rigid-body modes are least controllable 
from both locations as could be expected from the 


large numerical magnitude of the inertia matrix I s . 
Mode 1 is easier to control from location 1 than 2. 
Modes greater than 4 are less likely to experience 
spillover (ref. 5) effects from location 1. The rigid- 
body modes are most observable. Mode 1 is least 
observable from location 2 and actually behaves nu- 
merically as an unobservable mode in Kalman filter 
computations. Generally, table All shows that the 
first three modes have better controllability and ob- 
servability properties from location 1. 

The presence of invariant zeros (ref. 13) in the 
system models can be tested by seeking those values 
of s such that 


Rank 


sl n “ A f 
C p 


b f 

0 


< n + 3 


(A18) 


for given values of 1 < n q < 10. No zeros exist for 
n q — 0 since in this case 

det [C F (sl 6 - A f ) _1 B f ] = det (ij 1 ) ^ 0 

(A19) 

For n q > 0, the finite invariant zeros can be obtained 
from data given in table AIII. 


TABLE All. RECIPROCAL CONDITION NUMBERS FOR COLLOCATED SENSOR/ ACTUATOR 

LOCATIONS 1 AND 2 


Reciprocal 

Mode i 

Condition 

number 

Rigid 

body 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

C(«).(l) .... 

2.4 x 1(T 7 

4.0 x 10" 3 

1.7 x 10" 3 

1.9 x 10~ 3 

2.5 x 10" 6 

1.3 x 1CT 4 

5.6 x 10~ 5 

8.5 x 10" 5 

1.6 x 10“ 4 

4.1 x 10“ 4 

2.8 x 10“ 4 

C(«). (2) • • • 

2.4 x ltr 7 

1.9 x 10~ 5 

1.5 x 10" 4 

i.9 x icr 4 

5.9 x 10" 5 

2.3 x 1(T 3 

2.2 x 10" 3 

1.1 x 10~ 3 

6.4 x 10~ 4 

1.1 x 10~ 3 

5.5 x 10" 4 

000.(1) 

0.56 

1.6 x 1(T 3 

1.3 X 1(T 3 

1.5 x 10" 3 

2.4 X 1(T 6 

1.3 x 10“ 4 

5.5 x 10“ 5 

8.4 x 10" 5 

1.6 x 10“ 4 

4.0 x 10" 4 

2.7 x 10” 4 

000.(2) • • • 

0.56 

7.9 X 1(T 6 

1.1 x 1(T 4 

1.6 x lO -4 

5.6 x 10~ 5 

2.3 x 10 -3 

2.2 x 10 -3 

1.1 x 10“ 3 

6.3 x 10 -4 

1.0 x 10 -3 

5.5 x 10" 4 


TABLE AIII. DATA FOR INVARIANT ZEROS 


n q 

Sensor/actuator location 1 

Sensor/ actuator location 2 

1 

-0.9 x 10 -4 ±i0.082 

-0.75 x l(r^ ±/0.75 

2 

-0.37 x 10~ 3 ± j'0.22 

-0.013 ±yi.3 

3 

-0.29 x 1(T 3 ± j'0.22 

—0.017 ± jl.7 

4 

-0.032 ±i3.2 

-0.032 ±i3.2 

5 

—0.045 ± j4.5 

-0.79 x 10 -4 ±/0.19 

6 

-0.056 ± j’5.6 

-0.11 x 10 -3 ±/0.24 

7 

-0.057 ± ^5.7 

-0.17 x 10“ 3 ±j'0.32 

8 

-0.064 ± j’6.6 

-0.065 ± j6.5 

9 

-0.042 ±jb.l 

— 0.065 ± /6.5 

10 

-0.074 ±j7.7 

-0.082 ± j'8.3 
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Given a particular n q selection, the corresponding 
set of system zeros is the collection of entries from 
row 1 to row n q in table III. For example, with n q — 3 
and location 1, the set of system zeros are 

( -0.9 x 10~ 4 ) ± J.082 
(-0.37 x 1CT 3 ) ± j'0.22 
(-0.29 x 10~ 3 ) ± i0.22 


For both locations, all zeros are (marginally) min- 
imum phase. Many zeros lie close to the open-loop 
eigenvalues of A# (eq. (A5)) giving, as in table All, 
an indication that these modes will exhibit weak 
controllability /observability properties. Even though 
the zeros from location 2 for n q < 4 have higher fre- 
quencies (therefore allowing a higher controller band- 
width before interfering with performance), loca- 
tion 1 is chosen for the LQG/LTR design studies due 
to better controllability /observability properties. 
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TABLE I. ANTENNA PARAMETERS 


Weight, lb 

Inertia about axes through center of mass: 

I xx, lb-ft-sec 2 

Iyy, lb-ft-sec 2 

Izz, lb-ft-sec 2 

Structural mode frequencies, rad/sec . . 

Ixy = lyz = Ixz = 0 


10015.64 

4.222 x 10 6 

4.239 x 10 6 

3.233 x 10 6 

0 . 75 , 1 . 35 , 1 . 7 , 3 . 18 , 4 . 53 , 
5 . 59 , 5 . 78 , 6 . 84 , 7 . 4 , 8.78 


TABLE II. COMPENSATOR GAINS FOR ATTITUDE FEEDBACK CASE 


Filter gain matrix Hi 2 x 3 


4.47 

X 

10 “* 

- 9.32 

X 

10 -9 

4.36 

X 

10“ 8 

2.12 

X 

10-8 

4.47 

X 

10 -1 

7:72 

X 

10“ 8 

2.44 

X 

1Q-S 

1.48 

X 

10“ 7 

4.47 

X 

10“ 1 

1.00 

X 

10“ 1 

1.33 

X 

10“ 9 

7.60 

X 

10“ 9 

1.33 

X 

10“ 9 

1.00 

X 

10 “ J 

- 7.96 

X 

10“ 9 

7.60 

X 

10“ 9 

7.97 

X 

10“ 9 

1.00 

X 

10“ 1 

1.39 

X 

10“ 5 

1.42 

X 

10“ 5 

- 3.60 

X 

10“ 9 

7.92 

X 

10“ 6 

- 1.41 

X 

10“ 5 

- 1.94 

X 

10“ 9 

5.12 

X 

10“ 6 

1.43 

X 

10“ 6 

- 3.34 

X 

10“ 11 

2.59 

X 

10“ 6 

- 5.16 

X 

10“ 6 

1.19 

X 

10“ 10 

3.29 

X 

10“ 6 

5.69 

X 

10“ 7 

8.92 

X 

10“ 11 

1.65 

X 

10“ 6 

- 3.31 

X 

10“ 6 

1.26 

X 

10 -tO 


Control gain matrix G3X12 
10 5 

- 9.85 X 10“ 7 
7.61 X 10 2 
7.52 X 10 5 
1.23 x KT 5 
- 3.32 x 10 2 
- 2.63 
9.55 

2.57 x 10“ 8 
2.88 x 10“ 8 
1.14 x 10 2 
1.98 x 10 2 


- 2.73 x 10“ 7 
10 5 

4.5 x 10“ 7 
2.26 x 10“ 5 
7.61 x 10 5 
3.64 x 10“ 6 
- 1.91 x 10“ 10 
2.49 x 10“ 10 
- 9.31 x 10 1 
- 1.95 x 10 2 
- 6.81 x KT 9 
1.413 x 10“ 9 


- 7.61 x 10 2 
- 6.2 x 10“ 7 
10 5 

- 8.6 x 10 3 
- 7.44 x 10“ 6 
3.35 x 10 5 
- 4.12 x 10 2 
- 3.83 x 10 2 
3.67 x 10“ 9 
1.00 x 10“ 9 
3.60 x 10 1 
2.32 x 10 1 
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TABLE III. CLOSED-LOOP EIGENVALUES FOR ATTITUDE FEEDBACK 


Real part 

Imaginary part 

-8.54 x 1(T 3 

8.05 X lO -2 

-8.54 x 1(T 3 

-8.05 x lO -2 

-7.56 x 1(T 2 

1.25 x 10 _1 

-7.56 x 1(T 2 

-1.25 x 10 _1 

-7.60 x 10~ 2 

1.25 x HT 1 

-7.60 x 10“ 2 

-1.25 x nr 1 

-2.24 x 10“ 1 

2.24 x 10 -1 

-2.24 x 10“ 1 

-2.24 x lO -1 

-2.33 x 10 _1 

2.15 x 10 _1 

-2.33 X lO -1 

-2.15 X lO -1 

-2.38 x 10“ 1 

2.11 x 10 _1 

-2.38 X 10“ 1 

-2.11 x 10 _1 

-7.47 X 10 -3 

7.47 x 10 -1 

-7.47 x 10 -3 

-7.47 x 10" 1 

-1.35 x 10 -2 

1.35 x 10° 

-1.35 x 10“ 2 

-1.35 x 10° 

-3.08 x 10“ 1 

1.37 x 10° 

-3.08 x 10" 1 

-1.37 x 10° 

-1.02 x 10° 

1.27 x 10° 

-1.02 x 10° 

-1.27 x 10° 

-1.70 x 10“ 2 

1.70 x 10° 

-1.70 x 10 -2 

-1.70 x 10° 

-4.03 x 10- 1 

1.74 x 10° 

-4.03 x 10- 1 

-1.74 x 10° 

-3.18 x 10 -2 

3.18 x 10° 

-3.18 x 10- 2 

-3.18 x 10° 

—4.42 x 10 -2 

4.53 x 10° 

-4.42 x lO -2 

-4.53 x 10° 

-5.58 x 10~ 2 

5.59 x 10° 

-5.58 x 10“ 2 

-5.59 x 10° 

-5.73 x 10“ 2 

5.78 x 10° 

-5.73 x 10- 2 

-5.78 x 10° 

-6.69 x lO -2 

6.84 x 10° 

-6.69 x 10 -2 

-6.84 x 10° 

-6.39 X 10 -2 

7.40 x 10° 

-6.39 X 10" 2 

-7.40 x 10° 

-8.33 X 10“ 2 

8.78 x 10° 

-8.33 x 10- 2 

-8.78 X 10° 



TABLE IV. COMPENSATOR GAINS FOR ATTITUDE AND RATE FEEDBACK CASE 


Filter gain matrix Hi2x6 


4.17 X ur 1 

1.31 X KT 11 
1.67 X 10“ u 
9.09 X 10“ 2 
3.11 X HT 12 

2.32 x 10“ 12 
2.51 x 10“ n 
3.85 x 10“ n 

-1.67 x 10“ n 
8.07 x 10“ u 
6.40 x 10" 10 
2.83 x 10“ 10 


-2.87 x 10“ 12 

4.17 x 10 -1 
-9.32 x 10“ 12 
-2.16 X 10“ 13 

9.09 x 10“ 2 
-7.95 x 10“ 13 
2.47 X 10“ n 
-2.51 x 10“ n 
-7.10 x 10“ 10 
-3.32 x 10“ 10 
-1.85 x 10“ n 
9.03 x 10“ n 


1.51 x 10“ 12 
6.79 x 10“ 12 

4.17 x 10 _1 
1.84 x 10“ 13 
2.06 x 10“ 12 

9.09 x 10“ 2 
-3.32 x 10“ 9 
-1.83 x 10“ 9 
-3.15 x 10“ 11 
1.12 x 10 _1 ° 
8.03 x 10“ n 
1.19 x 10“ 10 


9.09 x 10“ 2 
3.92 x 10“ 12 
1.78 x 10“ 12 

4.17 x 10“ 2 
1.25 x 10“ 12 
2.63 x 10“ 14 

-7.10 x 10“ n 
3.37 x 10“ n 
-1.36 x 10“ 10 
7.91 x 10“ 12 
-5.92 x 10 -11 

2.18 x 10~ 9 


1.99 x 10“ 12 

9.09 x 10“ 2 
2.05 x 10~ 12 
7.07 x 10“ 13 
4.17 x 10“ 2 

6.09 x 10“ 13 
6.89 x 10“ n 
3.12 x 10“ 12 
6.88 x 10“ 11 

-1.58 x 10“ 9 
-7.70 x 10“ u 
-1.64 x 10“ 11 


-8.20 x 10“ 13 
9.36 x 10“ 13 

9.09 x 10“ 2 
-2.84 x 10“ 13 
4.52 X 10“ 13 
4.17 x 10“ 2 
5.39 x 10“ 10 
-2.91 x 10“ 9 
-8.74 x 10“ n 
-6.84 x 10“ 12 
-4.46 x 10“ u 
3.16 x 10“ 10 


Control gain matrix G ^ 2X 3 


10.00 x 10 3 

2.40 x 10“ 8 

-1.01 x 10 1 

5.05 x 10“ 9 

1.00 x 10 4 

-1.16 x 10“ 8 

1.01 x 10 1 

1.12 x 10“ 8 

10.00 x 10 3 

2.84 X 10 5 

1.39 x 10“ 6 

-4.57 x 10 2 

1.35 X 10“ 6 

2.85 X 10 5 

7.11 x 10“ 8 

4.00 x 10 1 

3.82 X 10“ 7 

2.11 x 10 5 

-7.10 x 10° 

-7.49 x 10“ 4 

-1.52 x 10° 

9.64 x 10“ 1 

2.40 x 10“ 14 

-7.87 x 10 1 

1.22 x 10“ 9 

-2.45 x 10° 

2.53 x 10“ u 

4.38 x 10“ 9 

-2.81 x 10 1 

-2.98 x 10“ 10 

3.13 x 10° 

7.35 x 10“ 10 

7.08 x 10“ 1 

3.50 x 10 1 

-2.53 x 10“ n 

5.19 x 10° 
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TABLE V. CLOSED-LOOP EIGENVALUES FOR ATTITUDE AND RATE FEEDBACK 


3.35 x 10~ 2 

3.35 x 10~ 2 

3.36 x 1(T 2 
3.36 x 10“ 2 
2.87 X 10 -2 
2.87 x HT 2 
2.32 X KT 1 

2.32 X 1CT 1 
2.30 X KT 1 
2.30 X 10 _1 
2.29 X 10 -1 
2.29 X 10“ 1 
7.47 X 10 -3 
7.47 X 10~ 3 
2.09 X HT 1 
2.09 X 10“ 1 
1.35 X 10“ 2 
1.35 x 10“ 2 
5.42 X 10 -2 
5.42 X 10" 2 
1.70 x 10" 2 
1.70 x 10“ 2 

8.33 x 10~ 2 
8.33 x 10 -2 
3.18 x 10 -2 
-3.18 x 10~ 2 
-4.51 x 10 -2 
-4.51 x 10 -2 
-5.59 x 10~ 2 
-5.59 x 10 -2 
-5.77 x 10 -2 
-5.77 x 10 -2 
-6.81 x 10~ 2 
-6.81 x 10" 2 
-7.19 x 10~ 2 
-7.19 X 10“ 2 
-8.78 X 10 -2 
-8.78 X 10 -2 


3.52 X 10“ 2 
-3.52 x 10 -2 

3.52 X 10~ 2 
-3.52 x 10~ 2 

4.47 x 10~ 2 
-4.47 x 10 -2 

2.14 X 10 _1 
-2.14 x 10 -1 

2.17 X 10 -1 
-2.17 x 10 -1 

2.18 x KT 1 
-2.18 x 10 -1 

7.47 x 10“ 1 
-7.47 X 10 -1 

7.53 X HT 1 
-7.53 x HT 1 

1.35 X 10° 
-1.35 X 10° 

1.35 x 10° 
-1.35 X 10° 

1.70 x 10° 
-1.70 x 10° 

1.70 x 10° 
-1.70 x 10° 

3.18 x 10° 
-3.18 x 10° 

4.53 x 10° 
-4.53 x 10° 

5.59 x 10° 
-5.59 x 10° 

5.78 x 10° 
-5.78 x 10° 

6.84 x 10° 
-6.84 x 10° 

7.41 x 10° 
-7.41 x 10° 

8.78 x 10° 
-8.78 x 10° 
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Feeds and feed panels 


1, 2 Actuator locations 



Figure 1. Hoop/column antenna concept. 







VSin 

Sanaa 



(a) Multiplicative uncertainty. 


Uncertainty 



(b) Additive uncertainty. 
Figure 3. Definition of uncertainty. 
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Singular values 










Figure 5. Block diagram of compensator/plant loop. 














11. Singular values of loop transfe 
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(b) Modified attitude and rate feedback loop. 

Figure 24. Modified block diagram for attitude and rate feedback. 
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